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Quantum input-output response analysis is a useful method for modeling the dynamics of complex 
quantum networks, such as those for communication or quantum control via cascade connections. 
Non-Markovian effects have not yet been studied in such networks. Here we extend the Markovian 
input-output network formalism developed in optical systems to non-Markovian cascaded networks 
which can be used, e.g., to analyze the input-output response of mesoscopic quantum networks. We 
use this formalism to explore the behavior of superconducting qubit networks, where we examine the 
effect of finite cavity bandwidths. We also discuss its application to open- and closed-loop control 
networks, and show how these networks create effective Hamiltonians for the controlled system. 

PACS numbers: 03.65.Yz,42.50.Lc,03.67.-a 



I. INTRODUCTION 

There has been tremendous progress in the last few 
years in experimental efforts to realize quantum net- 
works [l[ in various mesoscopic systems. These systems 
include photonic crystals 0], ion traps and super- 
conducting circuits [4|, which also advancesrelated fields 
such as quantum simulation (for recent reviews, see, e.g., 
Ref. (5|). The input-output formalism of Gardiner and 
coiiet |,g is a useful tool for analyzing such networks. 
In fact, using the input-output response to analyze or 
even modify the dynamics is a standard method in en- 
gineering called system synthesis. Up until now, sys- 
tem synthesis for quantum networks has only been stud- 
ied for Markovian systems. Quantum input-output the- 
ory itself has also mainly been limited to the Markovian 
regime, although it was developed for quantum systems 
about twenty years ago |g|. It was extended to the non- 
Markovian systems only quite recently [§[ . 

In the existing literature [H, 0] , quantum input-output 
theory is mainly applied to optical systems, in which the 
coupling between the system and its environment is weak 
and the correlation time (the "memory" of the environ- 
ment) is small compared with the characteristic time- 
scale of the system dynamics. Under the Markovian as- 
sumption, the quantum input-output formalism [6| was 
extended to cascaded systems Q, and has been used to 
study quantum feedforward and feedback networks [§- 
[llj . Markovian quantum input-output networks can 
be described using two alternative formulations: the 
Hudson-Parthasarity formalism in the Schrodinger pic- 
ture (12|; and the quantum transfer function formalism 
in the Heisenberg picture [H, [l4| • The general algebraic 
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structure of such systems has been well studied in the 
language of quantum Wiener and Poisson processes and 
quantum Ito rules [l2[. 

Although the Markovian assumption is reasonable 
when considering optical network components, environ- 
ments in mesoscopic solid-state systems can have corre- 
lations on much longer timescales [l5l-[Tcl|. Examples of 
this are the nuclear spin bath that couples to electron 
spins in quantum dots [l9j], and the 1/f noise that af- 
fects Josephson- Junction qubits [20]. It has also been 
suggested that the damping and decoherence of nano- 
mechanical resonators are due to coupling to a small 
number of two-level systems [2 if], which can be expected 
to induce significant non-Markovian dynamics. In addi- 
tion, any classical noise with a sufficiently narrow band 
generates non-Markovian evolution. There has been in- 
creasing interest in recent years in non-Markovian open 
quantum systems, and a number of analytical approaches 
have been devised to describe them. These include the 
projection-operator partitioning technique 123 , the non- 
Markovian quantum trajectory approach [231 ] . and very 
recently a non-Markovian input-output formalism Q . 

In this paper we extend the non-Markovian input- 
output theory to cascaded quantum networks, providing 
a recipe for obtaining non-Markovian input-output equa- 
tions for the description of any such network. Naturally, 
this formalism reduces to the standard input-output net- 
work formalism in the Markovian limit. Although the 
non-Markovian input-output relation has been derived in 
Ref. [§| , a quantum measurement is imposed on the out- 
put field in Ref. Q and the system dynamics is described 
by the quantum state diffusion equation. Such a for- 
malism cannot be extended to describe a non-Markovian 
cascade system, because quantum coherence in the out- 
put field is deteriorated by the measurement. In our 
formalism, without introducing measurements, the sys- 
tem dynamics is described by a non-Markovian quantum 



2 



stochastic differential equation and a perturbative mas- 
ter equation which can be naturally extended to a non- 
Markovian network. 

In developing our non-Markovian formalism, we will 
keep the weak-noise (weak-coupling) approximation, as 
this is the appropriate regime for implementing quan- 
tum technologies (for recent reviews, see, e.g., Ref. [24J). 
such as information processing and metrology. Here, we 
go beyond the Markovian approximation by allowing the 
coupling to the bath to have an arbitrary frequency de- 
pendence. This allows one to describe noise with any 
frequency profile, and should provide a good model for 
a wide range of non-Markovian environments. Interest- 
ingly, the resulting non-Markovian network formalism is 
exact for all couplings. However, in order to perform 
calculations for nonlinear systems, one must transform 
the Heisenberg equations of the input-output formalism 
to non-Markovian master equations, and this requires 
further approximations. Here we do this at the sim- 
plest level of approximation, by deriving the correspond- 
ing master equation to second-order in perturbation the- 
ory, using the standard Born approximation [22[. Nev- 
ertheless, more sophisticated techniques exist for obtain- 
ing non-Markovian master equations, and it would be 
an interesting avenue for future work to examine how 
these can be used to obtain master equations for non- 
Markovian cascaded networks. We note that for linear 
systems the Heisenberg equations of the input-output 
formalism can be used to obtain exact results. This is 
especially useful in some cases |25], when the second- 
order per turbative master equation fails to behave cor- 
rectly [26[ . We expect the network formalism we develop 
to be useful in describing a range of mesoscopic systems, 
such as coupled-cavity arrays in photonic crystals [27H281 ]. 
and nonlinear resonator and qubit networks in solid-state 
circuits [29-32]. The formalism can also be applied to 
quantum feedback control networks [33M39j | in solid-state 
systems [40^47}. 

This paper is organized as follows: in Sec.|Hl we briefly 
review the Markovian input-output formalism so that 
this can be easily compared to the non-Markovian case. 
In Sec. IIII1 we use an alternative method to derive the 
non-Markovian input-output relations in Ref. Q , and ob- 
tain the dynamical equation for the system such that it 
can be easily used for networks. We derive these here 
as a natural extension of the original Collett-Gardiner 
quantum input-output theory. In Sec. IIV1 we derive the 
input-output relations for more complex non-Markovian 
quantum cascade networks. We then apply these gen- 
eral results to two examples: a simple non-Markovian 
oscillator, and a network of two superconducting charge 
qubits interacting with Lorentz noises. We also apply 
the formalism to an open-loop and a closed-loop control 
networks, showing how these networks create effective 
Hamiltonians for the controlled system. The conclusions 
and prospects for future work are discussed in Sec. |Vl 



II. BRIEF REVIEW OF INPUT-OUTPUT 
THEORY OF MARKOVIAN SYSTEMS 

Here we summarize the standard Gardiner-Collet 
input-output formalism Q. The basic model is a quan- 
tum system interacting with a bath, where the bath con- 
sists of the modes of an electromagnetic field, or equiva- 
lently a continuum of harmonic oscillators. The Hamil- 
tonian for the system and bath is 



H — Hs + Hb + Hu 



(1) 



with 



Hr 



+ 00 



oj b*(uj) b(oj) dcu, 

> 

[k (uj) 6 t (u>) L — h.c] du>, (2) 



where w (oj) and b (oj) are the creation and annihilation 
operators of the bath mode with frequency oj, which sat- 
isfy 



[b (oj) , fe f (&)] =6(w-u>). 



(3) 



In the above, Hs is the free Hamiltonian of the system. 
The bath mode with frequency oj interacts with the sys- 
tem via the system operator L and the coupling strength 
K (oj). Hereafter we set H= 1. The total Hamiltonian H 
can be reexpressed in the interaction picture as 

H eS = exp (iH B t) (H s + H int ) exp (-iH B t) 



= H s + i 



-[-OC; 



[« (oj) e iut tf (oj) L - h.c] doj. (4) 



If the coupling strength is constant for all frequencies, so 
that 



K (oj) 



(5) 



then the dynamics of the system will become Markovian. 
This is the Markovian approximation. The Hamiltonian 
H e g is now given by 



H cS = H s +i^f \b\ n (t) L - Ltfe in (t) 



where 



hn (t) 



t b (oj) duj 



(6) 



(7) 



is the Fourier transform of the bath modes, and is the 
time-varying input field that is fed into the system [48[ 
(see Fig. [1]). In the Heisenberg picture, the system op- 
erator X (t) satisfies the following quantum stochastic 
differential equation (QSDE) 



\W [X,L] + [L\X]L} 



X = -i[X,H s 

+ ^{bin[L\X}+[X,L] bl}. 



(8) 
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FIG. 1: (Color online) Schematic diagram of the Markovian 
input-output system. 



If the input field is in a b- ln (i) vacuum state, and we 
trace it out, we can re-express the system dynamics in 
the Schrodinger picture as the following master equation 

p = -i [H S , P ]+j Uptf - l -tfhp - \ptfL^ . (9) 

Finally, if 6 ou t (t) is the field after it has interacted with 
the system, and is now out of it, then one has the relation 

b ou t(t)=b ia (t) + ^L(t). (10) 

This is the Markovian input-output relation. 



III. INPUT-OUTPUT THEORY OF 
NON-MARKOVIAN SYSTEMS 

A. General theory 

To derive the input-output relation for a general non- 
Markovian quantum system (see Fig. [2j , we rewrite the 
Hamiltonian H e g in Eq. (j4|) as 



H cS = H s + i 
= H s + i 



+ 00 



K(T-t)bl(r)dT 

bt(t)L-L% n (t) 



L - h.c. 



(11) 



where 



K(t) 



-foo 



exp (— iojt)n (w) duj (12) 



is the Fourier transform of the coupling strength k(uj). 
The input field that interacts directly with the system is 
now 



/+oo 
K(t - t) b in (r) dr 
-OO 

and satisfies the new commutation relation 

M*)»&L(*)1 =7 (*-*), 



where 



1 {t-t) 



+ 00 



[t — t) k (i — r) dr. 



(13) 



(14) 



(15) 



FIG. 2: (Color online) Schematic diagram of the non- 
Markovian input-output system. The input field is dispersed 
when it interacts with the system, and the modes with dif- 
ferent frequencies in the input field are coupled to the system 
with different coupling strengths. 



We can now proceed to derive the Heisenberg stochas- 
tic differential equations for the evolution of the system: 



X = -i[X,H s }+ /{ 7 (t-r)Lt(r)[A:(t),L(t)] 
Jo 

+7* (t-r) [Lt (t),X(t)]L(r)}dr 

+ {b in (t)[LKx]+[X,L]bl(t)}. (16) 

The non-Markovian input-output relation becomes 



(t) = 6 in (i) + K(t-r)L (r) dr. 



(17) 



Note so far that no further assumptions have been made, 
based on the Hamiltonian in Eq. (jlip . Thus our non- 
Markovian network formalism will give a more precise de- 
scription of non-Markovian systems with environment's 
coupling that have an arbitrary frequency-dependence. 
However, as noted in the introduction, for non-linear net- 
works we must often resort to the Schrodinger picture to 
perform calculations. 

To obtain the second-order perturbative master equa- 
tion, one averages over the vacuum input field b m , which 
we will take to be in the vacuum state, and uses the Born 
approximation. The perturbative master equation that 
corresponds to Eq. (IT^]) is 



-i [H S , p] + j o { 7 (* - t) [Lp (r) , 4 s (r - t) 
+7* {t-r) [L Hs (r-t),p(r)L^]}dr, 



where 



I/H s (t) = exp (iHst)Lexp [—iHst). 



(18) 
(19) 



We give the details of the derivations of Eqs. (TTBT) and 
(fT5]) . and output equation (fT?) in Appendix lAl 

Remark 1: The input-output relation, i.e., Eq . (fT7|) . 
coincides with Diosi's non-Markovian input-output 
equation (Eq. (10) in Ref. @). However, the dynamics 
of the input-output system in Ref. 0] is described by 
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FIG. 3: (Color online) Schematic diagram of the linear non- 
Markovian cavity. 



the quantum state diffusion equation after introducing 
a heterodyne measurement on the output field. Such a 
formalism cannot be extended to derive the dynamics 
of more complex non-Markovian networks, because 
quantum effects of the input-output system have been 
reduced by the quantum measurements. In our for- 
malism, without introducing measurement, the system 
dynamics is described by the non-Markovian quantum 
stochastic differential equation (|16l) and the perturbativc 
master equation (|18[) . Such a formalism can be more 
naturally extended to non-Markovian cascade networks. 

Remark 2: In the Markovian limit, where 

k (w) = \J^i/2it, we have n(t) — y^5(t), and 
7 (t — t) = j5 (t — t). It can be easily verified that 
Eq. (fT6|) reduces to the quantum stochastic differential 
equation ^ and the output equation (IT71) reduces to 
Eq. (TO]). 



B. Example: Single-mode cavity 

For a single-mode cavity coupled to an external input 
field (see Fig. [3]), the system Hamiltonian Hs and dissi- 
pation operator L are given by Hs = uj^a and L = a, 
where ujo, a (et^) are respectively the frequency and the 
annihilation (creation) operator of the cavity mode. The 
mode has an arbitrary non-Markovian coupling to the 
field modes outside the cavity. The quantum stochastic 
differential equation for a cavity operator X is then 



X 



-i [X, uj^a] + / { 7 (i - r) a) (r) [X (t) , a (t)] 
Jo 

+7* (t-r) [a* (t),X(t)]a(r)}dr 
+ {b- m [L^X]+[X,L]bl}. 



(20) 



Let us define the normalized position and momentum 



operators for the cavity as q = 



= a 



a 



/V2, p 



(— ia + m^) /v2- Collecting these into a single vector, 
x = (q,p) T , we may write the equations of motion for 



A x + Ai I 7 (t - t) x (t) dr 
Jo 

+A\ [ 7 * (t-r)x (r) dr + Bb m + B*b\ n , (21) 
Jo 



where 
A 



UJQ 
~0J 



-1 i 

-i -1 



(22) 



and B = {~l,i) T /V2. 

To solve the dynamical equation (j2"Tj) . we use the 
Laplace transform 



O(s) 



exp (-st)0 (t) dt 



(23) 



to transform the differential equation to an algebraic 
equation in the frequency domain. In this domain the 
solution becomes 



x(s) = 
where 



d(s) uj 
-w d(s) 
V2A(s) 



-K( S )b in (S)-K* (s)bl( S ) 

in (s) b- m (s) - in* (s) b[ n (s) 



(24) 



A(s) = d 2 (s)+u 2 , d(s) = s + j(s)/2, 



and k(s), x (s), b ln (s), bj n (s), 7 (s) are the Laplace trans- 
forms of k (t), x (t), bi n {t), bj n (t), 7 (t). We can also ob- 
tain the following input-output relation in the frequency 
domain 



"out 



iuj 



+ 7 (s) /2 + iuj 



(25) 



Note that this result is exact as far as the frequency 
dependence of the coupling to the bath is concerned. 
Because the system is linear we can obtain results 
without deriving a master equation. This input-output 
formula shows exactly how the coupling profile applies a 
low-pass filter to the input field to produce the output 
field. 

Remark 3: Let us consider a Markovian cavity with 
damping rate 7, then we have k (s) = w^y, 7 (s) = 7, 
from which we can obtain the traditional input-output 
relation for a lossy cavity from Eq. (|25|) (see Eq. (45) in 
Ref. Q) 



, , x s - 7/2 + iuj 

Oout («) = ; 777— h n [S) ■ 



s + 7/2 + iuj 



(26) 



IV. NON-MARKOVIAN QUANTUM 
NETWORKS 



A. General theory: quantum cascade systems 

To derive the input-output relation for complex non- 
Markovian networks, we should first study the dynamics 
of a system composed of two cascade-connected subsys- 
tems, also known as the "series product" of two subsys- 
tems [9]. The Hamiltonian of two cascaded subsystems, 
depicted in Fig. 2] can be expressed as 



System 1 



System 2 




FIG. 4: (Color online) Schematic diagram of the non- 
Markovian quantum cascade system. The output from the 
first subsystem is fed into the input of the second subsystem. 



H. 



off 



H x +i 
H 2 + i 



+ 00 



— oo 



«i (J - t) b\ M (r) dr 
n 2 (t - t) &J . (r) dr 



Li 
L 2 



h.c. 



where 



h.c. 



Jlr (t - t) = / k* (t - t) ki (i - r) dr. 



(32) 



(27) 



where ifi=i,2 and £j=i,2 are the free Hamiltonian and dis- 
sipation operator of the i-th subsystem; and Ki (t) is the 
corresponding coupling strength between the i-th subsys- 
tem and the i-th input field. If we omit the time delay for 
the quantum field transmitting between the two input- 
output components, then we have 



(t) = &i,out (*) = &i,in (*) + / «i (t - t) Li (t) dr 



Substituting Eq. fl5SJ) into Eq. ([27]). we have 
where 



(28) 



3,m 



i=l,2 



7l e 2 (r - t) L 2 L\ (t) - h.c. 



dr 



(29) 



is the interaction Hamiltonian between the two subsys- 
tems introduced by the transmitting field; the parameter 
J12 ( r — i s defined by 

/+oo 
k\ (t - t) k 2 (r - t) 6 (t - t) dr; 
-OO 

and 0[t) is the step function 

*<«>={ J; 1< 8: ( 3 °) 

The two equivalent non-Markovian input fields that inter- 
act directly with the two subsystems via the dissipation 
operators L\ and L 2 are defined as 

/+oo 
ki (t - t) b- m (t) dr, 
-OO 

and these satisfy the following commutation relation 

kta(t) ) Vj='Ylr{t-^, (31) 



The dynamics of the total system can then be expressed 
as the following quantum stochastic differential equation 

k = E f {lir(t-r)Lj(r)[X(t),L r (t)} 



l,r=l 



+lt(t-T) [4 (t),X(t)]L t (T)}dT 



It, 

-i[X,H s ] + ^2{h 



i=i 



L\,X 



[X,Lfbl n }, 
(33) 



where H$ = H\+ H 2 + H 12 . The input-output equation 
of the cascade system can be expressed as 



tout (*) = b- m (*) + / ki (t - t) Li (t) dr 



+ / «2 (* - t) L 2 (t) dr. 
Jo 



(34) 



Remark 4: In the Markovian limit, we have ki (t) = 
V^iS (*), 7f r (* - *) = 7/r (* - 1) = VTOr^ (*-*)> and 
&i,in = kin- Thus, the dynamical equation (|3T[) can be 
re-expressed as 



X = -i[X,H s ] + \{L^[X,L]+[L\X]L} 



b in [L\X]+[X,L] bl, 



(35) 



where 



Hq — Hi + H 2 + 



VTi72 



£ = y/jlLi + y^L 2 , 

and output equation (140[) can be rewritten as 

6out (*) = V (t) + (*) + V72"i2 (t) • (36) 

These equations coincide with those obtained for the 
Markovian series product systems in the literatures (see, 



6 



e.g., Ref. H). This concludes this remark. 



and Hi r is the field-mediated interaction Hamiltonian 



The dynamical equation (|33|) can be extended readily 
to n cascade-connected subsystems to obtain the follow- 
ing quantum stochastic differential equation 



* = E I \lir(t-r)Ll(r)[X(t),L r (t)} 

l,r=l J ° 

+Jt(t-r) [L\ (t),X(t)] U (r)}rfr 



7 b 

-i[X,H s ] + J2{ki 



1=1 



Lj,X 



+ [X,L l ] bh 



t \ 

Lin J ' 



H ir = i / % ( T ~ t) L r L\ (t) - h.c. dr. (39) 
Jo L J 

The function jf r (t — t) is defined by 

r+oo 

ll (*-*)=/ Kf 0" - t) Kr (t - t) (t - r) dr, 



where 8 (t) is the step function defined by Eq. (J30j) . The 
(37) output equation can be written as 



6out (i) = bin (i) + V / (t - r) L, (r) dr. (40) 
1=1 Jo 



where 7; r (t — f) is defined by Eq. (J32J ; 

n 

H s =Y. H i+Y. H ^ (38) 

(=1 i<r 

Transforming this into the Schrodinger picture, we can 
Hi = i_... tn is the free Hamiltonian of the Z-th subsystem; obtain the following second-order master equation 



J 




where 

£h s> ! ( f ) = ex P (iH s t)Li exp (~iH s t). 

B. Example: non-Markovian qubit networks in 
superconducting circuits 

As the first example, we apply our non-Markovian net- 
work formalism to superconducting circuits Here, as 
a simple example, we consider how to couple two distant 
single Cooper pair boxes (CPBs) by a microwave field. As 
shown in Fig. [SJ to suppress the decoherence effects, we 
embed two CPBs into two superconducting transmission 
line resonators (TLRs). When we average over the de- 
grees of freedom of the TLRs, the interactions between 
the CPBs and the input fields become non-Markovian. 
This can be understood by noting that the TLRs work 
as microwave cavities, and these act as low-pass filters. 
The white-noise input fields are filtered by the TLRs and 
changed into non-Markovian Lorentz noises, and these in 
turn interact with the CPBs. Thus, the qubit network 
considered here is a typical non-Markovian quantum net- 
work. 

The Hamiltonian of the j th CPB (j = 1,2) can be 
represented as 

ffcPBj = LE C (nj - n gj f - Ej ($ xj ) cos fa, (41) 



where <j>j denotes the phase drop across the j CPB, with 
rij = —id/ (d(j>j) as its conjugate operator. The operator 
rij represents the number of Cooper pairs on the island 
electrode. The scalar n g j is the reduced charge number 
on the control gate in units of Cooper pairs. This is given 
by n gj — ~CgVgj /2e, where C g and V g j are the gate ca- 
pacitance and gate voltage of the j th CPB. The scalar 
Ec = e 2 /2 (C g + 2Cj) is the single-electron charging en- 
ergy of a single CPB and Cj is the capacitance of a single 
Josephson junction. The Josephson energy Ej ($ X j) of 
the j th DC superconducting quantum interference device 
(SQUID) can be calculated by 

Ej($ xj ) = 2E j cos (ir^j , (42) 

where Ej represents the Josephson energy of a single 
Josephson junction; ^> x j is the external flux piercing the 
SQUID loop of the j th CPB; and $o is the flux quantum. 
For simplicity, we assume that Ec and Ej are the same 
for each Josephson junction in the two CPBs. 

Near the charge-degenerate point with n g j =0.5, the 
two lowest-energy levels of the j th CPB are close to each 
other and far separated from higher-energy levels. Thus, 
we can approximately consider a single CPB as a two- 
level system. In the qubit basis, the Hamiltonian -Hcpb.j 
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FIG. 5: (Color online) Schematic diagrams of (a) two cascade- 
connected CPB-TLR (Cooper-pair box and transmission line 
resonator) input-output systems and (b) the equivalent su- 
perconducting circuits. 



can be diagonalized as 



ffcPB, 



1„M 



(43) 



1/2 



E 



= E., 



where u q j 

E C j=4E c (l-2n gjJ . 

As shown in Fig. [Sj the j th CPB is capacitively coupled 
to the j th TLR. The Hamiltonian of the j th coupled CPB- 
TLR system can be written as 



Hi 



9j 



cosa,er 



00 



sin a 



where ay is the annihilation operator of the single- 



mode electromagnetic field in the j TLR; g = 
— e (Cg/Cs) Kms,j i s the coupling strength between the 
resonator and the qubit; Cs = C g + 2Cj is the total ca- 
pacitance of the SCB; V° ms = ^/u; c /2CV is the root mean 
square (rms) of the voltage across the LC circuit with C r 
representing the capacitance of the resonator; and aj is 
defined by 



oy = arctan [Ej (4> X j) / Ec (1 — 2n 



effective Hamiltonian of the j th coupled CPB-TLR sys- 
tem 



Hj = -Tp-crp' 1 + uJcjd^cij + gj + aj<r 



.00 



(44) 



where uj qj = Ej ((f> x j)- 

The cavity mode in the j th TLR is coupled to a trans- 
mitting field in an auxiliary transmission line between 
two TLRs. The total Hamiltonian of the j th SCB-TLR 
system and the input field b nitj can be represented by 



fftotj = +u) cj a]a j +g j (a] a 



00 



(45) 



where jj is determined by the coupling between the cav- 
ity mode and the input field. To be concentrated on the 
dynamics of the qubit, we eliminate the degrees of free- 
dom of the cavity mode. From Eq. (|4"3")) . it can be shown 
that 



lj 



lj b h 



igjo^. (46) 



Let us now introduce the following weak coupling as- 
sumption 

u c j, lj > gj, 

so we can omit the last term in Eq. (|46j) when we consider 
the dynamics of the cavity mode. Thus, we can now solve 
Eq. (06]) 



aj (t) 



Hj I exp[ 
lo 



(iujcj + 7j/2) (t - r)]6j,i n (r) dr. 

(47) 

(<£> •) and Substituting Eq. (|47|) into the Hamiltonian H to t.j in 



Eq. (|45|) , we can obtain an effective Hamiltonian to repre- 
sent the coupling between the j th qubit and the effective 
input field 



ft] 00 _ 00 £. . 



(48) 



where 



(*) = 



'ni-.i i ' ' — / l 9j 
lo 



/^ e -(i^+n/m-r) b . n . (t) dT _ (49) 



If we additionally add an ac gate voltage V g j = 
Voj cos {ojgjt) on the gate of the j th CPB, where Voj and 
ujgj are the amplitude and frequency of the gate voltage, 
we can obtain the following effective Hamiltonian in the 
rotating frame 



H e nj = ^J)+i($ n ^-afh 



(50) 



gj ■ 



under the condition that C g VojEc /2e <C A q j = 
Letting the SCBs work at the charge-degenerate point, Ej — us q j. By comparing Eq. (|4T))) and Eq. 
such that n g j = 1/2 (j = 1,2), and introducing the ro- we can see that the j CPB is just directly cou- 
tating wave approximation, we can obtain the following pled to the effective non-Markovian field bi n j with 



8 



Kj (t) — igjy/^fjexp [— (ui c j + jj/2) t}. Additionally we 
can see that the total system we consider here is just a 
cascade-connected two-qubit system mediated by a non- 
Markovian field. If &i„j is a white noise, it can be easily 
verified that the spectrum of 6; n j (t) is of Lorentz type. 
In fact, it can be calculated in the frequency domain that 



(u) = Kj H b in j (uj) 

jj/2 + i(u- uj cj ) 



M, (51) 



where &i n j (uj), 6; n ,j (u)), K (uj) are the Fourier transform 
of 6in,j (t), &i n .j (t), k (t). From Eq. (fT5")) . it can be shown 
that 



7j M = l«j M I 2 
Note that 
knJ (t),^j (t) 



7 2/4+( W - Wcj ) 



7j (* - *) 



J (* *)'yj (w) 



5 2 exp 



-7i|f-*| /2] 



We can see that 6i n j (t) is a Lorentz-type noise. 

To simplify our discussions, let us assume that the 
two qubits have the same system parameters, i.e., A g = 



^2, 9 



9i 



g 2 , and 7 



7i 



72 . From 



Eq. (|41[) . we can obtain the master equation of the non- 
Markovian two-qubit system we consider here 



if J z + a (t) a™*™ + a* (t) a^a^,p 



+ {I3(t-r) [J_p (r) , J+] + h.c.} dr, (52) 
Jo 



where J a =z,± = + o~a' is the collective two-qubit 
operator; and a (t) and j3 (t) can be calculated by 



r ( 2 ) 



:(t) = 



ig 2 exp(iA q t) ig 2 exp (— 7 t/2) 



7/2 + iAq 



7/2 + iA q 



(i{t) = g 2 exp [-"ft/2 + iA q t] 



(53) 



Additionally, if A g = Ej — ui q = 0, we can rewrite 
Eq. ()52|) by introducing the first Markovian approxima- 
tion, i.e., to replace p(r) in the last integral term in 
Eq. JUJ) by p(t), as 



p = -i 
where 



a(t) <T W a { + ) +h.c.,p 



a(t) = iT (<) =z— (1 

7 ■ 



T(t)V[J_]p, (54) 



,- 7 t/2 



and the superoperator P [J_] p is defined by 
L> [J_] p = J-p,J+ - \ j+J-P - \pJ+J- 



It can be verified that both a (t) and T (t) decrease when 
7 decreases. This means that both the coherent inter- 
action between the two qubits and the damping induced 
by the transmitting field decrease with the increase of 
the correlation time of the non-Markovian Lorentz noises 
hn,j (t) scaled by I/7. 

In Fig. [S] we show the evolution of the concurrence of 
the two qubits. The concurrence is defined by 

C (p) = max{Ai — A 2 — A 3 — A4, 0} . 

where p is the system density matrix given by Eq. (|54[) ; 
Aj's are the square roots of the eigenvalues, in a decreas- 
ing order, of the matrix p; and p* is the complex conju- 
gate of p. From Fig. [6] we see that the damping rate of 
the concurrence decreases as the bath coupling strength g 
decreases, but increases when the correlation time of the 
environment (r cnv = 1/7) increases. This means that 
in non-Markovian environments the two-qubit entangle- 
ment is preserved longer than in a Markovian environ- 
ment, which agrees with our intuition. 

The two-qubit dynamics given by Eq. (|52p can be 
extended to multi-qubit networks to study many-body 
physical phenomena, such as quantum entanglement and 
correlations. 



Example: feedforward and feedback control 
network 



In this subsection, we consider using networks for 
quantum control [50l455l |. A simple control network is 
composed of two cascade-connected subsystems. One is 
called the controller, and the other the plant. The pur- 
pose of the network is to improve the performance of the 
plant by connecting it to the controller. Whether such a 
strategy will work depends in general upon the properties 
of each. Such a configuration is useful if for some physi- 
cal reason the controller can be tailored in ways that the 
plant cannot . 

It is common to divide control networks into two 
classes. The first one is an open-loop or "feedforward" 
scenario in which the input field is first fed into the con- 
troller to obtain a control signal and then the control sig- 
nal is fed into the quantum plant to control the dynamics 
of the plant (see Fig. [7]). In this method, the control sig- 
nal contains no prior information of the system dynam- 
ics. The second one is a feedback control network |3(J - l49| . 
This is also called closed-loop control, in which the input 
field is first fed into the plant to extract the information 
we need. The information-bearing output field is then fed 
into the quantum controller to obtain an output control 
signal which is fed back to change the dynamics of the 
plant. Such a control system may be viewed as a three- 
partite quantum cascade system: plant-controller-plant. 
Since the signal fed into the controller contains real-time 
information about the state of the plant, we can use it to 
adjust the behavior of the controller. 
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FIG. 6: (color online) Evolution of the concurrence (a) for 
various non-Markovian correlation rates 7 (the blue curve, 
red triangle curve, and the green plus-sign curve correspond to 
7 = 0.2/t, 0.5/t, 1/t) and (b) for various coupling strengths 
g (the blue curve, red triangle curve, and the green plus-sign 
curve correspond to g = 0.1/r, 0.2/t, 0.3/r). Here r = 10 
ns. The decay of the concurrence speeds up when increasing 
the correlation rate 7 of the non-Markovian noises and the 
qubit-environment coupling strength g. 



By examining Eq. (|39p one can see that the quantum 
controller introduces Hamiltonian terms like 



H B = -1 



j(t-r) u c (t) cLt 



L + h.c. }, (55) 



where u c (t) = u c (X Ci \ (t) , ■ ■ ■ , X c , n (t)) is a function of 
the system variables X c 1 (£),••■ , X Cj „ (t) of the quan- 
tum controller and L is the system operator of the plant. 
The signal u c (t) is the control signal that is input to the 
plant. The idea is to design the controller dynamics so 
that the evolutions of X Ct i (t) , ■ ■ ■ , X c , n (t) generate the 
appropriate control signal, u c (t). The primary difference 
between feedforward (open-loop) and feedback (closed- 
loop) controls is that in the former the controller vari- 
ables X c ^i, ■ ■ ■ ,X c>n can be looked as extrogenous vari- 



Input J 



Controller Plant 
Control signal 




I Output 

1."' P~ 



(a) Quantum feedforward control 



Plant 



Input 



Control 



Controller 
(b) Quantum feedback control 



Output 



FIG. 7: (Color online) Schematic diagrams of (a) quantum 
feedforward control network and (b) quantum feedback con- 
trol network. Here, H c (H s ) and L c (L 3 ) are the Hamiltonian 
and dissipation operator for the controller (plant). 



ables that do not depend on the dynamics of the plant. In 
the latter, X Ct i, ■ ■ ■ , X C) „ are functions of the endogeous 
variables of the plant, and can be obtained in terms of 
the plant variables by solving the dynamical equation of 
the controller. These two control methods have differ- 
ent advantages, and the best performance might only be 
achieved by combing them. 

For a concrete example, let us choose the plant to be 
a single-mode cavity so that 



H s = io a o)a, L s = a, 



(56) 



and the controller to be a fully-controllable two-qubit 
system. The Hamiltonian and dissipation operator of 
the controller are 



fJ>d,cr+, 



L r = (T_ 



(57) 



Here A q is the detuning between the transition frequency 
of the qubit and the frequency of the extrogenous driving 
field and fid is a classical extrogenous control parameter. 

We now show what happens if we give the controller 
a fast damping rate, so that it adjusts very quickly to 
change in the plant. In this case we can average out the 
degrees of freedom of the controller (see the derivations 
in Appendix IB")) . For the case of feedforward control, the 
reduced Hamiltonian of the cavity can be written as 



Hi, 



w a aJa + [u* (t) a + u c (t) a'] 



(58) 
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Controller 



Cavity 




(a) Modifying the dynamics of the 
cavity bv feedforward control 



Cavity' 



Input 



Control 



Output 




Controller 

(b) Modifying the dynamics of the 
cavity by feedback control 



FIG. 8: (Color online) Modifying the dynamics of a cavity 
using a qubit inside the controller by (a) feedforward control 
and (b) feedback control. 



where u c (t) is a classical control parameter depending on 
the state of the qubit but does not depend on the dynam- 
ics of the controlled cavity mode a. This is simply the 
Hamiltonian of a driven harmonic oscillator. As a com- 
parison, for feedback control the reduced Hamiltonian of 
the cavity becomes 



H 



2,eff 



= uj a a)a + 



7* (t - t) a f (r) dr 



; (i) + h. 



(59) 

Non-classical optical effects such as squeezing induced by 
the memory term in Eq. (|59[) might well be observable in 
this case [56| ■ We see that feedback control can generate 
a class of evolutions that are impossible with an open- 
loop connection. 



V. CONCLUSIONS 

In summary, we have extended quantum input-output 
theory to arbitrary non-Markovian networks of systems 
connected via continuous-wave fields. We have derived 
the Heisenberg picture quantum stochastic differential 
equation for the systems in the network, the corre- 
sponding perturbative master equations, and all the 
input-output relations. We have applied this general 
formalism to a model of two superconducting charge 
qubits intracting via a cascade connection. We showed 



that this system was non-Markovian because the cavities 
with which the qubits connect to each other act as filters 
for the quantum noise. For this system we analyzed the 
dynamics of the entanglement between the qubits, and 
showed that it was affected by the non-Markovian nature 
of the network. We also use our model to analyze the 
difference between feedforward and feedback networks 
in which the controller has a fast response time. It is 
clear from our analysis that non-Markovian effects can 
have a significant effect on the behavior of mesoscopic 
quantum networks, and the analysis of these effects may 
be important for future quantum devices. 
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Appendix A: Derivations of the non-Markovian 
dynamical and output equations 

From the system Hamiltonian ([1]), we can obtain the 
Heisenberg equation of an arbitrary system operator X 

X = -i[X,H sys }+ J du{K(u)tf (u,t)[X,L] 

-«*(«) [X,tf]b(u,t)}, (Al) 
and the equation of the bath operator b (ui) 

b(u,t) = -iujb(uj,t) + k(lu) L. (A2) 
We can solve Eq. (|A2I) and obtain 

b (u, t) = e~ iut b (u) + k{uj) f e'^'^L (r) dr, (A3) 



where b(uj) = &(w,0) is the initial condition of b(u,t). 
Similarly, 



b(oj,t) = e- luj( - t - tl) b{w,t 1 )- K {u) I e'^-^L (r) dr, 

(A4) 
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where t\ > t. The input and output fields &i n (t) and 
bout (t) are defined as the Fourier transform of b (to) and 
b{u>,t\) respectively 



M*) = 

bout (t) = 



1 



/27T 
1 



-Hoc 



V2tt ,/___ 
From Eqs. (|A"3)) and (|A4]) . we have 



b (w) e^"* di, 





bout (t) = bin (t) + / ft (i — r) £ (r) dr. (A5) 
Jo 

Let ii — » i, we can obtain the output equation (|17|) . 
Furthermore, using the identities 



ft (t — t) &i n (t) dr 



— oo 



re (w) e^'b (w) dw, 



and 



7 (t-t) = 



+ 00 



k (w) ft* (w) e-M*-*) dw, 



we can obtain Eq. (fT75|) by substituting Eq. (|A3I) into 
Eq. (EH . 

To derive the master equation (fT8|). we first change 
into the interaction picture, in which the effective Hamil- 
tonian H c s can be rewritten as 



Hi, e g = i |&J n (t) £h s (t) - L Hs (t) bin (t) 



(A6) 



where Lh s (*) is given in Eq. (IT9l) . The density operator 
Pi, tot satisfies the following Liouville equation 



Pi,tot — —i [H^eS (t) , pi, tot] ■ 
Integrating the two sides of Eq. (|A7|) . we have 



(A7) 



Pi,tot (*) = -i [Hi, eS (t) , pi, e a (t)] dr. (A8) 



Substituting Eq. (|A8I) into Eq. (|A7[> . we can obtain 



Pi, tot 



[ffi, of f (t) , [i?i, ff (t) , pi,tot (r)]] dr. (A9) 



Tracing over the degrees of freedom of the input field, we 
can obtain the dynamical equation of the system density 
operator pj = tr B pi,tot 



Pi = I tr B {[-ffi.cff (t) , [-ffi, e ff (t) , pi,tot (r)]]} dr. 
Jo 



Let us then introduce the Born approximation and as- 
sume that the input field stays in the vacuum state, we 
have 

Pi,tot(i) = Pi{t)® |0) BB (0|. (All) 
Notice that it can be shown that 

(b- m (t)bt(t))=S(t-i), 

(bl (t) b in (<)) = (4 (t) fet n (<)) = (b in (t) & in (<)) = 0, 

(A12) 

where (•) is defined by (i?) = (0|b|0) _b- Substituting 
Eqs. (|A"o) . (|A11I) . and (|A12[) into Eq. (|A10|) . we can verify 
that 



Pi 



{ 7 (t - t) [L Hs (t) pi (r) , tf Hg (r)] + h.c.} . 



(A13) 

We can derive Eq. (fT8|) by transforming Eq. (|A13[) back 
into the Schrodinger picture. 



Appendix B: Derivations of the effective 
Hamiltonians for feedforward and feedback control 

For the case of feedforward control, the dynamics of 
the feedforward control system shown in Fig.JS^a) can be 
represented by 

<r_ = -iA q a- + i (p d + 9qbb) <r z , 

fr z = 2i (p d + g qb b) a- - 2i (p d + g qb b) a + , 

lb 

b = — -b-iuj b b + ig qb (j-+^/j b ~b- m , 

Kut = b in + y/jh'b, 

a = - + iw a ) a + VT^w, (Bl) 

where b is the annihilation operator of the cavity mode 
directly interacting with the qubit; -f a , lb are the decay 
rates of the cavity modes a, b; and g qb is the coupling 
strength between the qubit and the cavity mode b. Here 
we omit the decay of the qubit. By averaging out the 
input noise, we have 

b = ig qb / exp [- (iu b + j b /2) (t - r)] cr_ (r) dr. (B2) 
Jo 

Under the condition that pd 3> 9qb{b) where (b) is the 
average of the cavity mode b, we can rewrite Eq. (|Blj) by 
substituting Eq. (JEJ) into Eq. ([Bljl 



-iA q a- + ipdCFz, 



-{iu, b + lb /2)(t~T) 



a- (r) dr, 



(A10) 



bin + isjlbdqb I ' 
>0 

iVWb~9 qb [ e-^+^/ 2 )('- r V_(r)dr 
Jo 

- {la/2 + iuj a ) a. 



(B3) 



12 



Thus the dynamics of the cavity mode a is dominated by 
the effective Hamiltonian Hi, e g given in Eq. (|58[). where 

u c = -yjlalb9qb I exp [- (iui b + 7 6 /2) (f - r)] <r_ (r) dr. 
Jo 

Under the semiclassical approximation, er_ (t) can be re- 
placed by its average. With this simplification, u c (t) can 



be seen as a classical extrogenous control parameter be- 
cause the dynamics of the qubit determined by er_ and 
a z does not depend on the cavity mode a. With similar 
discussions, we can also obtain the effective Hamiltonian 
H 2 , c s induced by feedback control given in Eq. (l59|) . 
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